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Abstract 



In an exploratory study, using the tadpole improved clover Wilson quark action on 
small, coarse and anisotropic lattices, the irn scattering length in the 1 = 2 channel 
is calculated within quenched approximation. A new method is proposed which 
enables us to make chiral extrapolation of our lattice results without calculating 
the decay constant on the lattice. Finite volume and finite lattice spacing errors are 
analyzed and the results are extrapolated towards the infinite volume and continuum 
limit. Comparisons of our lattice results with the new experiment and the results 
from Chiral Perturbation Theory are made. Good agreements are found. 
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1 Introduction 



It has become clear that anisotropic, coarse lattices and improved lattice ac- 
tions are ideal candidates for lattice QCD calculations on small computers 
[1-4]. They are particularly advantageous for heavy objects like the glueballs, 
one meson states with nonzero spatial momenta and multi-meson states with 
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or without spatial momenta. The gauge action employed is the tadpole im- 
proved gluonic action on anisotropic lattices: 
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where Poi and Py represents the usual temporal and spatial plaquette vari- 
able, respectively. Rij and R i0 designates the 2x1 spatial and temporal Wilson 
loops, where, in order to eliminate the spurious states, we have restricted the 
coupling of fields in the temporal direction to be within one lattice spacing. 
The parameter u s , which is taken to be the fourth root of the average spa- 
tial plaquette value, is the tadpole improvement parameter determined self- 
consistently from the simulation. With this tadpole improvement factor, the 
renormalization of the anisotropy (or aspect ratio) £ = a s /a t will be small. 
Therefore, we will not differentiate the bare aspect ratio and the renormal- 
ized one in this work. Using this action, glueball and hadron spectra have 
been studied within quenched approximation [3-9]. In this paper, we would 
like to present our results on the pion-pion scattering lengths within quenched 
approximation. Main results of this paper has already been reported in [10]. 
Some details of the calculation are presented in this paper. 

Lattice calculations of pion scattering lengths have been performed by various 
authors using symmetric lattices without the improvement [11-13]. It turned 
out that, using the symmetric lattices and Wilson action, large lattices have 
to be simulated which require substantial amount of computing resources. It 
becomes even more challenging if the chiral, infinite volume and continuum 
limits are to be studied [13]. In this exploratory study, we would like to show 
that, such a calculation is feasible using relatively small lattices, typically of 
the size 8 3 40, with the tadpole improved anisotropic lattice actions. Attempts 
are also made to study the chiral, infinite volume and continuum limits ex- 
trapolation in a more systematic fashion. The final extrapolated result on 
the pion-pion scattering length in the 1 = 2 channel is compared with the 
new experimental result from E865 collaboration [14] and Chiral Perturbation 
Theory [15-18] with encouraging results. Of course, more accurate lattice cal- 
culations of the pion scattering length will inevitably require lattices larger 
than the ones studied in this paper. 

The fermion action used in this calculation is the tadpole improved clover 
Wilson action on anisotropic lattices [19,20]: 
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where various coefficients in the fermion matrix M. are given by: 

rji = v/(2u s ) , riQ = i/2 , a = \/{2k) - \/{\K max ) , 

p t = c 5w (l + 0/(4^) , p s = c sw /(2ut) . (3) 

Among the parameters which appear in the fermion matrix, csw is the coef- 
ficient of the clover term and v is the so-called bare velocity of light, which 
has to be tuned non-perturbatively using the single pion dispersion relations 
[20]. Tuning the clover coefficients non-perturbatively is a difficult procedure. 
In this work, only the tadpole improved tree-level value is used instead. 

In the fermion matrix (2), the bare quark mass dependence is singled out into 
the parameter a and the matrix A remains unchanged if the bare quark mass 
is varied. Therefore, the shifted structure of the matrix M. can be utilized to 
solve for quark propagators at various values of valance quark mass itlq (or 
equivalently k) at the cost of solving only the lightest valance quark mass 
value at k — « max , using the so-called Multi-mass Minimal Residual (M 3 R 
for short) algorithm [21-23]. This is particularly advantageous in a quenched 
calculation since one needs the results at various quark mass values to perform 
the chiral extrapolation. 

This paper is organized in the following manner. In Section 2, the method to 
calculate im scattering length is reviewed. In Section 3, some simulation details 
are described. In Section 4, our results of the scattering lengths obtained on 
lattices of various sizes and lattice spacings are extrapolated towards the chiral 
limit. A new quantity is proposed which is easier to measure on the lattice 
and has a simpler chiral behavior than the scattering length itself. Finite size 
effects are studied and two schemes of extrapolating to the infinite volume 
limit are studied. Finally, our lattice results are extrapolated towards the 
continuum limit. Comparisons with the experimental value and values from 
chiral perturbation theory at various orders are also discussed. In Section 5, 
we conclude with some general remarks. 



2 Formulation to extract the scattering length 



In order to calculate hadron scattering lengths on the lattice, or the scattering 
phase shifts in general, one uses Liischer's formula which relates the exact 
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energy level of two hadron states in a finite box to the scattering phase shift 
in the continuum. In the case of two pions which scatter at zero relative three 
momentum, this formula then relates the exact two pion energy E^ 1 ) in a finite 
box of size L and isospin / channel to the corresponding scattering length a 
in the continuum. This formula reads [24]: 
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where c\ = —2.837297, c 2 = 6.375183 are numerical constants. In this paper, 

(2) 

this formula will be utilized to calculate the pion-pion scattering length a in 
the 1 = 2 channel, which then requires the determination of the corresponding 
energy shift SE^J = — 2m 1T in that channel. 

One issue that one has to keep in mind is that, in a quenched calculation, 
Liischer's formula (4) has to be modified dramatically as discussed in [25]. Due 
to the anomalous contributions in quenched chiral perturbation theory from ?/, 
the energy shifts are contaminated by terms that might be of order L° = 1 and 
L~ 2 in the 1 = channel. If these contributions were substantial, extraction 
of the corresponding scattering length from the energy shift might become 
problematic. In the 1 = 2 channel, the situation is better. The quenched 
chiral corrections modify the energy shift by terms that are at least of order 
L~ 3 with usually small coefficients. Therefore, doing quenched calculation in 
the I = 2 channel is safer. 



To measure the pion mass and to extract the energy shift 5E<$ of two pions 
with zero relative momentum, we construct the correlation functions from the 
corresponding operators in the 1 = 2 channel. We have used the operators 
proposed in Ref. [12]. First, the local operators which create a single pion 
with appropriate isospin values are: 



7r + (x, t) = -d(x,t)7 5 w(x,t) , 7r (x,t) = w(x,t)7 5 (i(x,t) , 

A*, *) = t)75u(x, t) - J(x, t)7 5 rf(x, t)} , (5) 

where w(x, t), d(x, t), -u(x, t) and J(x, t) are basic local quark fields correspond- 
ing to u and d quarks respectively. The operators which correspond to zero 
momentum single pions are: 

4(0 = 75/2 XXfr*) ' ^ 
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where the flavor of pions a = +, — , and L 3 is the three volume of the lattice. 
Zero momentum one pion correlation function can then be formed as: 

c«(t) =< K(t)K(o) > . (7) 

From the large t behavior of this correlation function, the pion mass m n is 
obtained. Similarly, we use the two pion operators in the 1 = 2 channel defined 
by: 

I = 2 (t)=4(t)n+(t + l) , (8) 

to construct the two pion correlation functions: 

C T = 2 (t) =< Oif(t)Oif(0) > . (9) 

From the large t behavior of C^ 2 (t), one could infer the exact two pion state 
energy in the finite box. Numerically, it is more advantageous to construct the 
ratio of the correlation functions defined above: 

n I= \t) = ci=\t)/{c^t)c^t)) . (io) 

This ratio thus exhibits the following asymptotic behavior for large t: 

nit) t> ^ 1 e- SE ^ , (11) 

with 8E<® = EW - 2m w is the energy shift in this channel which directly 
enters Liischer's formula (4). It is argued [26,12,25] that this formula can only 
be utilized for small SE^t values in a quenched calculation for large enough 
L. In this case, one could equally use the linear fitting function: 

H(t) T>> ^ >M 1 - 6E$t , (12) 
to determine the energy shift SE^J. 

Two pion correlation function, or equivalently, the ratio H(t) constructed 
above can be transformed into products of quark propagators using Wick's 
theorem. In order to avoid the complicated Fierz re-arrangement terms, we 
used the creation operators at time slices that are differ by one lattice spacing 
as suggested in [12]. It can be shown that the 1 = 2 two pion correlation 
function is given by two contributions which are termed Direct and Cross con- 
tributions [11,12]. These contributions are represented by diagrams as shown 
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Fig. 1. Diagrams representing the contributions to the two-pion correlation function 
in the 1=2 channel. Diagrams (A) and (D) correspond to the Direct contribution 
while (B) and (C) belong to the Cross contribution. 

in Fig. 1, where the solid and dashed lines represent the u and d quark prop- 
agators respectively.QThe two pion correlation function in the 1 = channel 
is, however, more complicated which involves vacuum diagrams that require 
to compute the quark propagators for wall sources placed at every time-slice, 
a procedure which is more time-consuming than the 1 = 2 channel. In the 
1 = 2 channel, it turns out that the quark propagators have to be solved only 
twice with two fixed wall sources placed at t — and t — 1. 



3 Simulation details 

Simulations are performed on several PC's and two workstations. Configura- 
tions are generated using the pure gauge action (1) for 4 3 40, 6 3 40 and 8 3 40 
lattices with the gauge coupling (3 = 1.7, 2.2, 2.4 and 2.6. The spatial lattice 
spacing a s is roughly between 0.18fm and 0.39fm while the physical size of the 
lattice ranges from 0.7fm to 3.2fm. For each set of parameters, several hun- 
dred decorrelated gauge field configurations are used to measure the fermionic 
quantities. Statistical errors are all analyzed using the usual jack-knife method. 
Basic information of these configurations is summarized in Table 1. 

Quark propagators are measured using the Multi-mass Minimal Residue algo- 
rithm for 5 different values of bare quark mass. Periodic boundary condition is 
applied to all three spatial directions while in the temporal direction, Dirichlet 
boundary condition is utilized. In this calculation, it is advantageous to use 

1 In this work, the u and d quark are assumed to be degenerate in mass. 
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Table 1 

Simulation parameters for lattices studied in this work. Input aspect ratio parameter 
£ is fixed to be 5 for all lattices being studied. The approximate spatial lattice 
spacing a s in physical units as obtained from [4,9] is also indicated. Also listed are 
the maximum value of the hopping parameter K max - 



(3 Lattice u A s 


a s (fm) 


No. confs 


''max 


V 


1.7 4 3 x 40 0.295 


0.39 


700 


0.0570 


0.90 


1.7 6 3 x 40 0.295 


0.39 


288 


0.0585 


0.90 


1.7 8 3 x 40 0.295 


0.39 


240 


0.0585 


0.90 


2.2 4 3 x 40 0.378 


0.27 


700 


0.0590 


0.95 


2.2 6 3 x 40 0.378 


0.27 


256 


0.0600 


0.95 


2.2 8 3 x 40 0.378 


0.27 


224 


0.0600 


0.95 


2.4 4 3 x 40 0.409 


0.22 


700 


0.0620 


0.87 


2.4 6 3 x 40 0.409 


0.22 


256 


0.0605 


0.92 


2.4 8 3 x 40 0.409 


0.22 


224 


0.0605 


0.93 


2.6 4 3 x 40 0.438 


0.19 


600 


0.0630 


0.80 


2.6 6 3 x 40 0.438 


0.19 


256 


0.0620 


0.88 


2.6 8 3 x 40 0.438 


0.19 


240 


0.0605 


0.93 



the wall sources which greatly enhance the signal [11-13]. Values of the maxi- 
mum hopping parameter K max , which corresponds to the lowest valance quark 
mass, are also listed in Table 1. Typically, a few hundred Minimal Residual 
iterations are needed to obtain the solution vector for a given source vector. 
On small lattices, in particular those with a low value of ft, the hopping pa- 
rameter has to be kept relatively far away from the critical kappa value in 
order to avoid the appearance of exceptional configurations. 

The single pseudo-scalar and vector meson, which will be referred to as pion 
and rho, correlation functions at zero spatial momentum and three lowest 
lattice momenta, namely (100), (110 and (111) are constructed from the cor- 
responding quark propagators. Conventional effective mass functions are com- 
puted from which the single meson energy levels with zero and non-zero mo- 
mentum are extracted from the effective mass plateaus. The fitting range of 
the plateau is determined automatically by requiring the minimum of the x 2 
per degree of freedom of the fit. Using the anisotropic lattices, we are able 
to obtain decent effective mass plateaus from the single pion and rho corre- 
lation functions. In Fig. 2, the effective mass plateaus for the pion and rho 
are shown for one of our parameter set, namely 8 3 40 lattices with (3 = 1.7. 
Points with error bars are the effective mass values from the pion (lower half) 
and rho (upper half) correlation functions. For simplicity, only data for one 
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Fig. 2. The single pseudo-scalar (lower) and vector (upper) meson effective mass 
plot obtained from the meson correlation functions on 8 3 40 lattices with (3 = 1.7. 
Only the results for one value of k is shown. The bare velocity of light is taken to 
be v = 0.90. The open triangles, open squares, open pentagons and open hexagons 
are data for the zero and the lowest three non-zero spatial lattice momenta ((000), 
(100), (110) and (111)), respectively. Grey shaded horizontal bars represent the 
fitted values, the errors and the fitting ranges of the corresponding energy levels, as 
explained in the text. 

particular value of k is shown. The grey shaded horizontal regions represent 
the final fitted value of the energy. The starting and ending positions of these 
shaded regions correspond to the fitting ranges. The height and thickness of 
each of these shaded region designates the fitted value and the corresponding 
error respectively. It is seen that single meson energy values are obtained with 
good accuracy, especially for the mass. The fitting quality for other parameter 
sets is similar. The accuracy of the mass values are so good that we could ne- 
glect the errors of these mass values when we perform the error analysis of the 
scattering length where the errors from the energy shifts SE^ are dominating. 

Having obtained the single meson energy levels for zero and non-zero spatial 
momentum, the dispersion relation of the meson is also at our disposal. The 
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Fig. 3. The dispersion relation for the pseudo-scalar meson obtained on 8 3 40 lattices 
with (3 = 1.7. The bare velocity of light is taken to be v = 0.90. It is seen that the 
measured dispersion relations agree with the free dispersion relation at small lattice 
momenta with a slope of 1.06(3). 

parameter u, also known as the bare velocity of light, that enters the fermion 
matrix (2) is determined non-perturbatively using the single pion dispersion 
relations as described in Ref. [20]. In Fig. 3, we show the single pion dispersion 
relation obtained at the optimal value of v — 0.90 for (3 = 1.7, 8 3 40 lattices. 
Points with error bars are the single pion energy levels obtained by fitting the 
pion effective mass plateaus. All data which correspond to 5 different values of 
k are plotted. The straight lines are the linear fits to these data. The optimal 
value of v is such that the slope of these lines are consistent with unity. The 
optimal value of v which corresponds to each parameter set of the simulation 
is also tabulated in Table 1 

Two pion correlation functions and the ratio TZ(t) are constructed from prod- 
ucts of suitable quark propagators according to Wick's theorem. Both the 
Direct and the Cross contributions are included. For the ratio lZ(t), we ob- 
tain good signal for all our data sets. In Fig. 4 we show the ratio 1Z(t) for 
all five k values as a function of the temporal separation t, together with the 
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Fig. 4. The ratio R(t) as a function of t obtained from the two pion correlation 
functions at zero lattice momentum for 5 different values of the bare quark mass. 
The lattice size is 8 3 40 and the gauge coupling f3 = 2.2. The optimal bare velocity of 
light is taken to be v = 0.95. The straight lines are linear fits according to Eq. (12), 

(2) 

from which the energy shifts 5E = — 2m n are extracted for all 5 values of k. 

corresponding linear fit (12) for (3 = 2.2 on 8 3 40 lattices. The fitting range 
cannot start at a value of t that is too small, particularly in a quenched calcu- 
lation, as discussed in [25]. The starting and ending positions of the straight 
lines indicate the appropriate fitting range which minimizes the x 2 per degree 
of freedom. Here, especially on lattices with larger sizes and better signal, 
we have tried to enlarge the fitting range as much as possible. But for small 
lattices, e.g. 4 3 40 lattices at f3 — 2.6, signals only last a few time slices. 

Finally, the extracted values of the single meson energy and the energy shifts 
are all listed in Table 2,Table 3 and Table 4 for 4 3 40, 6 3 40 and 8 3 40 lattices, 
respectively. 

After obtaining the energy shifts SE$, these values are substituted into 

(2) 

Liischer's formula to solve for the scattering length % for all five k values. 
This is done for lattices of all sizes being simulated and for all values of (3. The 
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Table 2 

Extracted values of the single pseudo-scalar, vector meson mass and the energy shift 
6E^I in the isospin 1 = 2 channel obtained from the simulation on 4 3 40 lattices, all 

(2) 2 

measured in aV 1 . Also listed are the values of the quantity F = a ° mp (see Sec. 4 



for the discussion) obtained by substituting the value of SE^J into Eq. (4). 





P = 


1.7, 4 3 40 lattices 




(3 = 2.2, 4 3 40 lattices 


K 


m w 


m p 


5E (2) 

(10-4) 


—F 


K 


m n 


m p 


5E (2) 

(10-4) 


—F 


.057 


.4129(6) 


.512(2) 


99(3) 


1.29(3) 


.059 


.275(1) 


.334(2) 


200(3) 


1.04(1) 


.056 


.4827(6) 


.570(2) 


89(3) 


1.42(3) 


.058 


.341(1) 


.392(2) 


175(3) 


1.24(1) 


.055 


.5503(5) 


.629(1) 


81(3) 


1.57(4) 


.057 


.406(1) 


.451(1) 


157(3) 


1.44(1) 


.054 


.6163(5) 


.689(1) 


75(3) 


1.72(4) 


.056 


.4700(9) 


.510(1) 


143(2) 


1.66(2) 


.053 


.6813(5) 


.7496(8) 


69(2) 


1.88(5) 


.055 


.5330(8) 


.570(1) 


131(2) 


1.88(2) 




= 


2.4, 4 3 40 lattices 




(3 = 2.6, 4 3 40 lattices 


K 


m n 


m p 


6E<$ 

(10-4) 


—F 


K 


m w 


m p 


(10-4) 


—F 


.062 


.180(2) 


.216(3) 


256(19) 


0.58(3) 


.063 


.215(1) 


.246(2) 


143(15) 


0.45(4) 


.061 


.245(1) 


.276(3) 


218(17) 


0.78(4) 


.062 


.2701(9) 


.286(2) 


135(13) 


0.56(4) 


.060 


.308(1) 


.336(2) 


186(18) 


0.97(6) 


.061 


.3261(8) 


.342(2) 


125(12) 


0.73(5) 


.059 


.368(1) 


.395(2) 


140(10) 


1.03(5) 


.060 


.3821(8) 


.398(2) 


116(11) 


0.89(6) 


.058 


.428(1) 


.453(2) 


158(21) 


1.45(12) 


.059 


.4380(7) 


.454(1) 


107(10) 


1.06(7) 



results for the quantity F = a m \/m^ (see the next section for the reason of 
this choice) are also listed in Table 2 to Table 4. From these results, attempts 
are made to perform an extrapolation towards the chiral, infinite volume and 
zero lattice spacing limit. 



4 Extrapolation towards chiral, infinite volume and continuum lim- 
its 



Since the valance quark mass values being studied are far from the chiral limit, 
we first try to make the chiral extrapolation at finite volume and finite lattice 
spacing. This is facilitated by results for the scattering length at 5 different 
values of k. In the chiral limit, the nn scattering length in the 1 = 2 channel 
is given by the current algebra result due to S. Weinberg [15]: 

a (2) - - — ^ (U) 
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Table 3 



Same as Fig. 2 but for 6 3 40 lattices. 



(5=1.7, 6 3 40 lattices 


(3 = 2.2, 6 3 40 lattices 


K 


77V 


m p 


5E {2) 

UJ - / 7T7T 

(10-4) 


—F 


K 


77V 


m p 


5E {2) 

UJ - / 7T7T 

(10-4) 


—F 


.0585 


.3029(8) 


.416(2) 


43(6) 


1.32(14) 


.060 


.202(1) 


.294(2) 


81(6) 


1.19(7) 


.0575 


.3803(7) 


.477(1) 


37(5) 


1.45(15) 


.059 


.2751(7) 


.348(2) 


66(5) 


1.33(8) 


.0565 


.4472(5) 


.537(1) 


32(4) 


1.62(16) 


.058 


.3440(8) 


.404(1) 


56(5) 


1.51(9) 


.0555 


.5159(4) 


.597(1) 


30(4) 


1.82(18) 


.057 


.4097(7) 


.461(1) 


49(5) 


1.72(12) 


.0545 


.5827(4) 


.657(1) 


28(4) 


2.03(20) 


.056 


.4737(6) 


.520(1) 


43(4) 


1.91(13) 


(3 = 2.4, 6 3 40 lattices 


(3 = 2.6, 6 3 40 lattices 


K 


m n 


777 p 


«5£ (2) 

7Y7Y 

(10-4) 


—F 


K 


77V 


777 p 


UJ - / 7Y7Y 

(10-4) 


—F 


.0605 


.205(1) 


.267(3) 


66(21) 


0.82(20) 


.0615 


.199(1) 


.232(2) 


130(10) 


1.09(6) 


.0595 


.275(1) 


.322(2) 


62(8) 


1.08(10) 


.0605 


.2626(9) 


.291(1) 


106(8) 


1.38(6) 


.0585 


.3402(9) 


.379(1) 


48(10) 


1.18(19) 


.0595 


.3241(6) 


.351(1) 


91(7) 


1.70(8) 


.0575 


.4030(8) 


.436(1) 


38(10) 


1.25(25) 


.0585 


.3840(5) 


.409(2) 


79(6) 


2.00(10) 


.0565 


.4610(6) 


.494(2) 


31(10) 


1.30(32) 


.0575 


.4439(6) 


.467(1) 


71(6) 


2.30(12) 



where a is the Tin scattering length ao in the 1 = 2 channel and ~ 93MeV 
is the pion decay constant. Chiral Perturbation Theory to one-loop order gives 
an expression which includes the next-to- leading order contributions [16]. Even 
the next-next-to-leading order corrections have been calculated within Chiral 
Perturbation Theory [17]. For a summary of the results in Chiral Perturbation 
Theory, see Ref. [18] and the references therein. Since the one- loop and two- 
loop numerical results on the pion-pion scattering length in the 1 = 2 chan- 
nel do not differ from the current algebra value substantially, we will mainly 
compare with Weinberg's value. Complication arises in the quenched approx- 
imation. In principle, the quenched scattering lengths becomes divergent in 
the chiral limit [25] . However, these divergent terms only become numerically 
important when the pion mass is close to zero. For the parameters used in our 
simulation, these terms seem to be numerically small in the 1 = 2 channel and 
we could not observe this divergence from our data. 

In most of the previous lattice calculations [11-13], both the mass and the 
decay constant of the pion were calculated on the lattice. Then, the lattice 
results of m n and f n were substituted into the current algebra result (13) 
to obtain a prediction of the scattering length. This is to be compared with 
the scattering length obtained from Luscher's formula. In these studies, some 
discrepancies between the lattice results and the chiral results was observed 
[13]. There were also discrepancies between the staggered fermion results and 



12 



Table 4 

Same as Fig. 2 but for 8 3 40 lattices. 



/3 = 1.7, 8 3 40 lattices 


P = 2.2, 8 3 40 lattices 


K 




m p 


5E {2) 

7Y7Y 

(10-4) 


—F 


K 


m n 


m p 


5E {2) 

UJ - / 7Y7Y 

(10-4) 


—F 


.0585 


.3005(7) 


.421(2) 


21(3) 


1.57(21) 


.060 


.2026(7) 
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P = 2.4, 8 3 40 lattices 


P = 2.6, 8 3 40 lattices 
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the Wilson fermion results. The major disadvantage of the this procedure is 
the following: First, it is difficult to make a direct chiral extrapolation towards 
the chiral limit; Second, the lattice results of the decay constant are usually 
much less accurate, both statistically and systematically, than the mass values. 
In fact, most of the discrepancies are due to inaccuracy of the decay constants, 
as realized already in [13]. In this paper, we propose to use another quantity 
Q which has a simpler chiral behavior to compare with the current algebra 



result (13). This quantity is simply F 
reads: 



rrip/m^, which in the chiral limit 



(2) 2 



1 m p 



-1.3638 



(14) 



where the final numerical value is obtained by substituting in the experimental 
values for m p ~ 770MeV and f % ~ 93MeV. Note that both m p and f n have 
finite values when approaching the chiral limit. Therefore, the above combi- 
nation will approach a finite value in the chiral limit. On the lattice, we could 

(2,) 

extract the scattering length a from Liischer's formula. More importantly, 
the mass of the pion and the rho meson can be obtained with good accuracy 

2 A similar dimensionful, instead of dimensionless, quantity a^/m^has been pro- 
posed in Ref.[13]. 
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on the lattice. So, the factor F can be calculated on the lattice with good pre- 
cision without the lattice calculation of f n . The error of the factor F obtained 

(2) 

on the lattice will mainly come from the error of the scattering length Oq , or 
equivalently, the energy shift 8E^} and nowhere else. Since we have calculated 
the factor F for 5 different values of valance quark mass, we could make a chi- 
ral extrapolation and extract the result in the chiral limit. Comparisons with 
Weinberg's result (14) and the experiment will offer us a cross check among 
different methods. 
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Fig. 5. Chiral extrapolation for the quantity F = m 2 p /m^ for our simulation 
results at (3 = 2.2 on 6 3 40 lattices. In the lower half of the plot, the pseudo-scalar 
mass squared are plotted as open squares as function of 1/(2k). The straight line 
represents the corresponding linear fit (15) for the data. Also shown in the lower 
half as open hexagons are the results for the factor F. In the upper half of the plot, 
the same quantity F is plotted (open squares) as a function of m%. Also plotted 
as open circles are the data points for m? p . The straight line represents the linear 
extrapolation towards the chiral limit = 0, where the extrapolated result is also 
depicted as a solid square. As a comparison, the corresponding experimental result 
[14] for this quantity is drawn as a shaded band. Weinberg's result and the result 
from chiral perturbation theory are also shown as a filled triangle and a filled circle 
at m n = 0, respectively. 

14 





In the upper half of Fig. 5, we show the chiral extrapolation of the quantity 
F = a^vn^/m^ as a function of the pseudo-scalar mass squared (m 2 ) for 
the simulation on 6 3 40 lattices at (5 — 2.2. The open squares represent data 
points for F at various values of m 2 . The straight line is the linear chiral 
extrapolation and the extrapolated result is also depicted as a solid square 
at m n = 0. Weinberg's result and the result from chiral perturbation theory 
are also shown as a filled triangle and a filled circle at m n = 0, respectively. 
The vector meson mass values squared are also shown as the open circles in 
the upper half of this plot. The pseudo-scalar meson mass squared (open 
squares) are plotted in the lower half of the figure as a function of l/(2/c), 
which linearly depends on the valance quark mass m q via: 

m « = A il-ir) , as) 



It is seen that meson mass squared depends on the valance quark mass linearly, 
as expected from chiral symmetry. It is also possible to dig out the critical 
value of the hopping parameter where the pseudo-scalar mass vanishes. The 
data points for the factor F are also shown in the lower half of the plot as open 
circles. The fitting quality for the pion, rho and the factor F are reasonable. 
The quality for the chiral extrapolation of other simulation points are similar. 
As is seen, the linear fit gives a reasonable modeling of the data. The divergent 
contributions from quenched chiral perturbation theory seem to be numerically 
small for the lattices being simulated. 

After the chiral extrapolation, we now turn to study the finite volume effects 
of the simulation. According to formula (4), the quantity F obtained from 
finite lattices differ from its infinite volume value by corrections of the form 
1/L 3 . However, it was argued in Ref. [26,12,25] that in a quenched calculation, 
the form of Luscher's formula is invalidated. Finite volume corrections will 
have a different dependence on the volume, e.g. a correction of the form 1/L 5 
instead of 1/L 6 , as predicted by formula (4). This would mean that the fac- 
tor F receives finite volume correction of the form 1/L 2 . In our simulation, 
however, we were unable to judge from our data which extrapolation is more 
convincing. We therefore perform our infinite volume extrapolation in both 
ways, calling them scheme I ( extrapolating according to 1/L 3 ) and scheme 
II ( extrapolating according to 1/L 2 ), respectively. In fact, extrapolation in 
these two different schemes yields compatible results within statistical errors. 
The fitting quality of Scheme II is somewhat, but not overwhelmingly, better 
than that of Scheme I. In Fig. 6 and Fig. 7, we have shown the infinite vol- 
ume extrapolation according to Scheme I and Scheme II, respectively, for the 
simulation points at /3 — 2.6, 2.4, 2.2 and 1.7. The extrapolated results are 
shown with open squares at L = oo, together with the corresponding errors. 
The straight lines represent the linear extrapolation in 1/L 3 or 1/L 2 . It is seen 
that, on physically small lattices, e.g. those with (3 = 2.6 and (3 = 2.4, the 
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0.02 
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Fig. 6. Infinite volume extrapolation in Scheme I for the quantity F = a m^/m^ 
for our simulation results at (3 = 2.6, 2.4, 2.2 and 1.7. The straight lines represent the 
corresponding linear extrapolation in 1/L 3 . The extrapolated result is also shown, 
together with its error. 

finite volume correction is much more significant than larger lattices. For the 
lattices with (3 = 1.7, the finite volume dependence of the result is rather weak, 
showing that even on 4 3 40 lattices with this lattice spacing, the finite volume 
correction is small. This is also reflected by the slopes of the linear fits. It is 
seen from the figures that the slopes of the fitting straight lines increase with 
(3, showing a larger finite volume correction at larger (3 values hence smaller 
lattices. 



Finally, we can make an extrapolation towards the continuum limit by elim- 
inating the finite lattice spacing errors. Since we have used the tadpole im- 
proved clover Wilson action, all physical quantities differ from their continuum 
counterparts by terms that are proportional to a s . The physical value of a s for 
each value of f3 can be found from Ref. [4,9], which is also included in Table 1. 
This extrapolation is shown in Fig. 8 where the results from the chiral and 
infinite volume extrapolation discussed above are indicated as data points in 
the plot for all 4 values of (3 that have been simulated. In the lower/upper 
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half of the plot, results in Scheme I/II are shown. The straight lines show the 
extrapolation towards the a s = limit and the extrapolated results are also 
shown as solid squares together with the experimental result from Ref. [14] 
which is shown as the shaded band. For comparison with chiral perturbation 
theory, Weinberg's result (14) and the results from Chiral perturbation theory 
are also shown as the filled triangles and filled circles at a s = 0, respectively. 
It is seen that our lattice calculation gives a compatible result for the quantity 
a< o* >m ' 2 p/ m n when compared with the experiment. The statistical errors for the 
final result is still somewhat large. This is mainly due to lack of results at 
smaller lattice spacings. However, the result is promising since the chiral, infi- 
nite volume and continxuum extrapolation for the pion scattering length have 
not been studied systematically using small lattices before. Another encour- 
aging sign is that all data points at different lattice spacing values are about 
the same which indicates that the 0(a s ) lattice effects are small, presumably 
due to the tadpole improvement of the action. To summarize, we obtain from 
the linear extrapolation the following result for the quantity F = a^m 2 Jm^. 
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Fig. 8. Continuum extrapolation for the quantity a mi/m^ obtained from our 
simulation results at (3 = 2.6, 2.4, 2.2 and 1.7. The results for both Scheme I and 
II are shown. The straight lines represent the linear extrapolation in a s /rQ. The 
extrapolated results are also shown, together with the experimental result from 
Ref. [14] indicated by the shaded band. For comparison, Weinberg's result (14) and 
the results from Chiral perturbation theory are also shown as the points at a s = 
respectively. 

-£ = -1.05(23) for Scheme I , 



rn 

,( 2 )™2 



a ° mp = -1.41(28) for Scheme II . 



(16) 



If we substitute in the mass of the mesons from the experiment: m p ~ 770MeV 

(2) 

and ~ 139MeV, we obtain the quantity % m^. 



a (2) m n = -0.0342(75) for Scheme I , 
ai 2) m w = -0.0459(91) for Scheme II . 



(17) 
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Weinberg's current algebra prediction (13) yields a value of a ( 'm n = —0.046. 
This quantity has been calculated in Chiral Perturbation Theory to one-loop 
order with the result: a^m^ = —0.042 [16] and recently to two-loop order 
[17,18]. The final result from Chiral Perturbation Theory gives: offm^ = 
—0.0444(10), where the error comes from theoretical uncertainties. On the 

(2) 

experimental side, a new result from E865 collaboration claims a m w = 
—0.036(9). It is encouraging to find out that our lattice results in both schemes 
are compatible with the experiment. Our result in Scheme II also agrees with 
the Chiral Perturbation Theory results very well while our result in Scheme I 
is barely within one standard deviation of the chiral results. 



5 Conclusions 



In this paper, we have calculated pion-pion scattering lengths in isospin 1 = 2 
channel using quenched lattice QCD. It is shown that such a calculation is fea- 
sible using relatively small, coarse and anisotropic lattices with limited com- 
puter resources like several personal computers and workstations. The calcu- 
lation is done using the tadpole improved clover Wilson action on anisotropic 
lattices. Simulations are performed on lattices with various sizes, ranging from 
0.8fm to about 3fm and with different value of lattice spacing. Quark propa- 
gators are measured with 5 different valance quark mass values. These enable 
us to explore the finite volume errors and the finite lattice spacing errors in a 
more systematic fashion. The infinite volume extrapolation is explored in two 
different schemes which yields compatible final results. The lattice result for 
the scattering length is extrapolated towards the chiral and continuum limit 
where a result consistent with the experiment and the Chiral Perturbation 
Theory is found. We believe that, using the method described in this ex- 
ploratory study, more reliable and accurate results on pion scattering lengths 
could be obtained with simulations on larger lattices and computers. 

Finally, our method for calculating the pion-pion scattering length discussed 
in this paper can easily be generalized to calculate the scattering lengths of 
other hadrons, or in other channels, e.g. 1 = channel, where extra care has 
to be taken due to enhanced terms coming from quenched chiral loops. The 
method can also be applied to calculate the scattering phase shift at non- 
zero spatial lattice momenta, where presumably lattices larger than the ones 
utilized in this exploratory study is needed. 
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